TOTALLY ASYMMETRIC ZERO-RANGE PROCESS 
IN THE RAREFACTION FAN 



PATRICIA GONCALVES 

ABSTRACT. We consider the one-dimensional totally asymmetric zero-range process starting 
from a step decreasing profile leading in the hydrodynamic limit to the rarefaction fan of the as- 
£N| . sociate hydrodynamic equation. Under that initial condition, we show that the weighted sum of 

joint probabilities for second class particles sharing the same site, is convergent and we compute 
its limit. We derive the Law of Large Numbers for the position of a second class particle initially 
at the origin under the initial state in which all positive sites are empty and all negative sites 
are occupied and also for a slight perturbation of the invariant state. 
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1. Introduction 



Interacting particle systems were introduced in the early seventies by Spitzer 1161 and in 
the hydrodynamic limit theory one obtains the partial differential equations governing the 
evolution of the conserved thermodynamical quantities of the system. Usually, one can take 
advantage of the known results on the hydrodynamic equation to obtain knowledge on the 
asymptotic behavior of the underlying particle system. Nevertheless, the study of the particle 
^ ! system has given answers on the qualitative behavior of the solutions of the corresponding 

■ hydrodynamic equation. In this paper we continue the work initiated in [ 5 ] - namely we study 

the asymptotic behavior of second class particles added to a particle system and relate them 
with the characteristics of the corresponding hydrodynamic equation. We are interested in 
the rarefaction fan setting and we prove that among the infinite characteristics emanating 
from the position of the second class particle, this particle chooses at random one of them. 
The randomness is given in terms of the weak solution of the hydrodynamic equation through 
some sort of renormalization function. 

In this paper we consider a well know interacting particle system: the one-dimensional 
totally asymmetric zero-range process (tazrp) evolving on Z. The dynamics of this process is 
the following: at each site there is a mean one exponential time clock, after which a particle 
^ ■ at x jumps to x + 1 at rate one independently from the number of particles at the destination 

site. After that jump the clock restarts. For this process, the equilibrium measures that 
are translation invariant are Geometric product measures of parameter with p e [0, oo), 
that we denote by fi p . Since we are restricted to the one-dimensional setting, we can couple 
this zero-range process with the totally asymmetric simple exclusion process (tasep) that was 
analyzed in |5 ]. The dynamics of the latter process is as follows: after a mean one exponential 
time, a particle at a site x jumps to x + 1 at rate one, only if the destination site is empty, 
otherwise it does not move. The invariant measures for this process are Bernoulli product 
measures of parameter a e [0, 1], that we denote by v a . By coupling both tazrp and tasep, we 
will be able to confirm the results that we prove independently for the tazrp, from the results 
in©. 

Since the work in Ifl4ll it is known that, the empirical measures associated to these pro- 
cesses, converge to a deterministic measure whose density is the unique entropy solution of 
a hyperbolic conservation law with concave flux. The hydrodynamic limit for these processes 
was set in two different ways: for both processes in 1 14] using the Entropy method, for a 
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general set of initial measures associated to a profile p and for the tasep in (8) using the 
Relative Entropy method, for a more restricted set of initial measures. We notice that it is 
not difficult to show the hydrodynamic limit for the tazrp invoking the same arguments as 
in (8] to derive its hydrodynamic limit. With last result in hands, we can also cover the case 
of initial measures of slowly varying parameter associated to flat profiles as well as the case 
of step decreasing profiles. 

The main goal of this paper is to analyze the asymptotic behavior of the tazrp starting 
from an initial configuration with particles with different degree of class. We mainly consider 
first or second class particles, but we define the interaction dynamics between particles for 
all m degree of class with m = l,..,oo. Initially each particle is labeled as a 771-th class 
particle and its movement depends on its label, in such a way that a m-th class particle 
sees particles with degree of class less (greater) than m as particles (holes). We can think 
of a hole as a particle with 00 degree of class. Now, we describe the interaction dynamics 
between particles of different degree of class sharing the same site. Suppose that at the 
site x there are particles of degree of class 1, ..,00. If the clock at x rings, the first class 
particle has priority to jump and moves to the right site at rate one, while the remaining 
particles keep the same position. So, if at a site, there is no particle with degree of class 
less than m, the 777-th class particle jumps to the right at rate one. The higher the degree 
of class of a particle, the lower is its priority to jump. On the other hand, for the tasep the 
interaction dynamics between particles with different degree of class is completely different. 
In the latter process, a m-th class particle can move forward to the right neighboring site if 
the destination site is occupied with a particle with degree of class greater than m, and in 
this case, they interchange positions. Nevertheless, differently from the tazrp, a m-th class 
particle can jump backwards if at the left neighboring site there is a particle with degree of 
class less than m, that attempts to jump to the right. So, in the tasep a m-th class particle 
interchanges with a particle with degree of class less than 777 and moves to the left with rate 
one, or interchanges with a particle with degree of class greater than m and moves to the 
right at rate one. 

Now, we recall some results that have been proved for the tasep, which we will generalize 
for the tazrp in a similar setting. Fix p, A e [0, 1] and let v Pt \ = v P &(-oofl] + ^a<5(o,oo)- In © it 
was shown, that starting the tasep from v Pi \ with p > A and adding a second class particle 
at the origin, it holds a Law of Large Numbers (L.L.N.) for the position of this particle and 
the limit IX is uniformly distributed on [1 — 2p, 1 — 2A]. Later, this result was extended to 
more general partially asymmetric transition rates of exclusion type in |4]. The speed of the 
second class particle in the partial asymmetric simple exclusion was also studied in Q] by 
analyzing the invariant measures of the multi-class process. It was proved in II 131 . that this 
L.L.N, holds in the strong sense; and for p = 1 and A = in this result was also derived 
by mapping the tasep to a last passage percolation model. As in [6] the second-class particle 
can be viewed as an interface between two random growing clusters, see |4] for details. 

For the tasep, the relation between the distribution function F u (u) and the entropy solu- 
tion p(t, u) of the associate hydrodynamic equation under initial condition Pq X (u) := pl{ u <o} + 
-^l{ji>o} is given by F u (u) = p ^ p _^ u ' ■ The purpose of this work is to see how far this picture 
can go in the tazrp. 

Let p Pj \ — /-t P <5(_oo,o] + /- < a<5(o,oo), Pi A € [0,oo]. Distribute initially the tazrp according to 
Moo,o- By adding a second class particle at the origin, a L.L.N, can be derived and the limit 
X is distributed according to F x (u) = 1 — 4>(p(l,u)), where 4>(p) — is the mean of the 
instantaneous current (see section ETcfl ) and p{t,u) is the unique entropy solution of the hy- 
drodynamic equation 02.11 * with initial condition p^' x . Be aware that, in this case the entropy 
solution does not define a density of a probability distribution function, so that, the relation 
between F x and p is performed by means of a proper function cf>. The idea of the proof of last 
result, consists in defining the current through a time-dependent bond and applying a cou- 
pling argument introduced in |5|. This argument seems robust enough to obtain the L.L.N. 
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for second class particles for attractive particle systems with hyperbolic hydrodynamic equa- 
tion with concave flux as in ( 12. 2D . On the other hand, considering the tazrp starting from (i Pt \ 
with p > A and adding infinite second class particles at the origin at the initial time, we are 
able to show that the sum of the joint probability of the speeds of these particles, converge 
and we compute its limit. For the tasep, the joint distribution of the speeds of second class 
particles was analyzed in Q]. 

As in the case of the tasep, we can also recover a L.L.N, for the position of the second class 
particle in the tazrp, in which the limit distribution function is given as F 2 (u) = p ~^^ , 

where pit, u) is the unique entropy solution of ( 12. ID with initial condition pft . For that pur- 
pose, we consider the tazrp under a slight perturbation of the invariant state p\. first we 
distribute particles according to and then we add a first class particle to each negative 
site and put a second class particle at the origin. 

Finally, we notice that by coupling the tazrp with the tasep and invoking Theorem 2.3 of 
0), we show that starting the tazrp from a configuration with a second class particle at the 
origin, a third class particle at the site 1, all negative sites occupied by infinite first class 
particles and all positive sites empty, the probability that the second class particle overtakes 
the third class particle is 2/3. It would be an interesting problem to derive this result for the 
tazrp without going through the coupling with the tasep. 

Here follows an outline of this paper. In the second section, we introduce the totally asym- 
metric zero-range process, we state its hydrodynamic limit, we compute the characteristics 
of hyperbolic conservation laws with concave flux and state the main results. On the third 
section, we compute the limit of the sum of crossing probabilities for the tazrp starting from 
fJ- P ,x (p > A) and with infinite second class particles at the origin. On the fourth section, 
starting the tazrp from ^oo,o, we prove the L.L.N, for the position of the second class particle 
and for the current of first class particles that cross over the second class particle. On the 
fifth section, starting the tazrp from a perturbation of the invariant state p, p (which is still 
associated to Pq' x ), by adding a second class particle at the origin, we derive its L.L.N. The 
sixth section, is devoted to the reproof of the results of the fourth section, under the coupling 
with the tasep and invoking the results in[5]. Here we also discuss the crossing probability 
of a second class particle to overtake a third class particle in the tazrp as in (4). 

2. Statement of results 

2.1. The dynamics. The one-dimensional tazrp is a continuous time Markov process £ t with 
state space N z . In this process, after a mean one exponential time, a particle at x jumps to 
x + 1 at rate 1, independently from the number of particles at the destination site. For 
a configuration £ and for x e Z, £(x) denotes the number of particles at that site, i.e. if 
= k, k = 0, .., oo, then there are k particles at the site x. Its infinitesimal generator is 
defined on local functions / : N z — > R as 

mo = J2 1 i^)>i}[f(^ x+1 )-f(0], 

where 

if z 7^ x, x + 1 
t(x)-l, itz = x 
+ + ifz = x + l 

We will also consider more general totally asymmetric jumps given by g(k), where k = 
and the rate function g : N — > K + is non-decreasing, g(Q) = 0, g(k) > for k > and with 
bounded variation, i.e. sup fc \g(k + 1) — g(k)\ < +oo. 

Now we introduce the dynamics of the tasep. Let ry t be a continuous time Markov process 
with space state {0, 1} Z . In this process, particles evolve on Z according to interacting random 
walks with an exclusion rule which prevents to have more than a particle per site. The 
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dynamics is as follows: after a mean one exponential time, a particle at x jumps to x + 1 at 
rate 1, if the destination site is empty otherwise it does not move and the clock restarts. For 
a configuration rj and for x G Z, tj(x) denotes the quantity of particles at x: rj(x) = 1, the 
site x is occupied otherwise it is empty Its infinitesimal generator is given on local functions 
/ : {0,1} Z -> K by n/(r?) - E* eZ ^)(l - f(s + l))[/(*f ,x+1 ) - /fa)] where rf'^ 1 is the 
configuration obtained from rj exchanging the variables at site x and x + 1. We notice that in 
both dynamics described above, particles only move in the one-dimensional lattice, so that, 
the total number of particles is a conserved quantity for both processes. 

Now, we recall the notion of attractiveness of the zero-range process. In N z there is a 
partial order between configurations as follows: C < C if Vx e Z ((x) < ({x). This partial 
order induces the corresponding stochastic order on the distributions of the process: if p and 
p are probability measures in N z , then there exists a coupling measure p in the product space 
N z x N z , whose marginals are p and p and supported in {((, C) : C < C}- It is well known 
that the fact of <?(•) begin non-decreasing implies that given £ and C such that C < C. then it 
is possible to construct a coupling of the zero-range process (Ct, Ct) starting from (£, £) such 
that Vi > Ct < Ct, see (2)> or Ifl2ll . Since we took <?(•) non-decreasing, we have in hand 
the attractiveness property for the processes we are considering. Later we will see that this 
property is crucial for our conclusions. 

Now we describe a set of equilibrium measures that we will recall in the sequel. For the 
tazrp, for each density p there exists an invariant measure p p , translation invariant and 
such that E Pp [£ (0)] = p, namely the Geometric product measure of parameter j±- : 

p p (Z:£(x) = k)= -J-,fce0,..,oo. 
VI + pJ 1 + p 

For the tasep, it is well known that for < a < 1, v a the Bernoulli product measure on 
{0, 1} Z with density a: 

v a {r] : i](x) = 1) = a, 

is an invariant measure, translation invariant and also such that E Va [?7(0)] = a. 
Let pp^\ be the product measure m p £(-oo,o] + p\5(o,oo), namely for k e 0, .., oo: 

\ (ji-Y-n-, itx<0 
{ (tTa) TTA' iix>0 
and denote by v p .\ the product measure ^ p <5(-oo,o] + ^a^o.oo) • 

2.2. Hydrodynamic limit. Now we state the hydrodynamic limit for the processes intro- 
duced above. Fixed £, let 7r"(£, du) be the empirical measure given by 

7r"(e,^) = ^^)5 f (rfu), 

where 5 U denotes the Dirac measure at u and let 7r"(£, du) = n n (^ t , du). 

Since the work of Rezakhanlou II 1411 . it is known that, starting the tazrp from a measure p n 
associated to a profile p and some additional hypotheses (for details see [ 14]), if 7Tq (£, du) con- 
verges to p (u)du in p n -probability, then 7r™ n (£, du) converges to p(t, u)du in p n S tn -probability, 
where S t is the semigroup corresponding to £ and p(t, u) is the unique entropy solution of the 
hyperbolic conservation law: 

d t p(t,u)+V<t>(p(t,u)) = 0, (2.1) 
where cf)(p) = For the process with jump rate <?(•) satisfying the assumptions above, the 
hydrodynamic limit also holds and in this case, 4>(p) = g(p) = E Pp [5(77(0))]. 

In 11141 the hydrodynamic limit was derived for the tasep and its hydrodynamic equation is 
the inviscid Burgers equation, i.e. equation Q2.1H with cf>(p) = p(\ — p). In [81 the hydrodynamic 
limit for the tasep was derived, via the Relative Entropy method and the same arguments can 
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be applied to the tazrp to derive its hydrodynamic limit under this method. As a consequence 
the local equilibrium convergence is readily obtained, see Corollary 6.1.3 of H10I . 

2.3. Characteristics. Now, we go towards a description of characteristics for hyperbolic 
conservations law as above, as in flU . Consider a partial differential equation given by: 

d t p(t,u) + VF(p(t,u)) = 0, (2.2) 

where the flux F(-) is a concave function and suppose it is differentiable. A characteristic 
is the trajectory of a point with constant density and denoting by v Pn (t, u) the position at 
time t of a point with density p = p(0, u), then p(v Po (t, u), t) = po = p(0, u). Taking the time 
derivative of last expression it follows that d s v pa (s, u) = F'(po). Integrating last expression 
from time to time t and noticing that w Po (0, it) = u, we get to v Po (t, u) = u + F'(p )t. So, in 
this case, characteristics are straight lines with slop F'(po). 
If the initial condition is a decreasing step function as 

Po X ( u ) = P l {u<o}(u) + Al {u>0 }(u), (2.3) 
with p > A, then the solution of i2.2i can be explicitly computed and it is given by 

( p, \£u<F'(p)t 
P {t,u) = \ A > tfu>F'{X)t . 
1^ ip(^j^j, otherwise 

The function ip can be easily computed and reads as tjj(v) = the inverse of F'. 

Now we derive these solutions for the processes under consideration. We notice that </>(p) 
corresponds to mean with respect to the invariant state of the process of the instantaneous 
current at a bond [x, x + 1], The instantaneous current at the bond [x, x + 1] is defined as the 
difference between the jump rate to the right neighboring site and the jump rate to the left 
neighboring site. In our cases, since jumps are totally asymmetric, it coincides with the jump 
rate to the right. For the tazrp, the instantaneous current at the bond [x, x + 1] is l{^(a;)>i}, 
which in turns gives 4>(p) = So, under the initial condition p p ' A , the solution at time t of 

( EO t is given by as above with F'{-) := <£'(•) and ii{u/t) := 

For the zero-range process with jump rate given by <?(•), the instantaneous current at 
the bond [x,x + 1] is g(^(x)) and the solution of the corresponding hydrodynamic equation 
under initial condition pfr x , is given on u € [(^'(p)) -1 *, (g'(A)) _1 t] by g'{j where g{p) = 
E Vp [3(77(0))] and outside of this interval is given as above. 

2.4. Discrepancy. Now we introduce the notion of discrepancy between two copies of a 
same attractive process. For that purpose, let r\. and £. be two copies of a process such that 
iJo(x) = £0(2) for all x 7^ and for example at the site zero take 770 (0) = £o(0) + 1. So, at 
time zero there is only one discrepancy between 770 and £ that is located at the origin. By 
the attractiveness property and the conservation of the number of particles, at each time t 
there is still a unique discrepancy between r) t and £ t . This discrepancy is called second class 
particle. As mentioned in the introduction, the dynamics of second class particles is very 
different in zero-range and exclusion type dynamics. On the tazrp, when first and second 
class particles share the same site, the first class particles have priority to jump, a second 
class particle jumps to the right neighboring site if and only if there is no first class particles 
at the departure point. On the other hand, for the exclusion type dynamics, the second class 
particle jumps to the right neighboring site if it is empty, but if a first class particle attempts 
to jump to the site occupied by the second class particle, they exchange positions. So, in the 
exclusion type dynamics a second class particle can jump backwards, contrarily to the tazrp 
in which it only moves to the right. Throughout the paper we denote by X 2 (t) the position 
of the second class particle on the tazrp at time t. Naturally, that copying two tazrp, we 
do not have necessarily one discrepancy between them, e.g. if we take in the copy above 
?7o(0) = £o(0) + k, with k e N, then there are k second class particles at the origin at time 0. 
The dynamics between second class particles sharing the same site is fully described below 
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(see the proof of Theorem |2.1D , but the idea is simple: second class particles at the same site 
are labeled from the bottom to the top, the first to jump to the right neighboring site is the 
one at the bottom and it only jumps if there is no first class particles at the departure point. 
For the tazrp we show that the sum of joint distributions of these discrepancies converge 
and we identify its limit. The joint distribution of the speeds of second class particles for 
the tasep was studied in [1]. Here we provide the first result on joint distributions of these 
speeds for the tazrp. 



2.5. Main results. 

Theorem 2.1. Consider the tazrp starting from p P! \ with < A < p < oo and suppose that 
at the initial time there are infinite second class particles and no first class particles at the 
origin. At the initial time label the second class particles, from the bottom to the top and for 
j € 2, .., oo, let Xj(t) be the position at time t of the j — th second class particle initially at the 
origin. Then 

Jim E > «*)[7^7 - WTW}) = + ' + - 

j—Q v ' v ' 

where pit, u) is the entropy solution of (2.1} with initial condition Pq X as in (2.3} . 

In the same spirit as in |5], we establish the L.L.N, for a single second class particle 
initially at the origin for the tazrp starting from Po^q. This is the content of next theorem. 

Theorem 2.2. Consider the tazrp starting from poo.o and at the initial time put a second class 
particle at the origin and remove the first class particles there. Denote by X 2 it) the position at 
time t of the second class particle initially at the origin. Then 

lim = X, in distribution 

t 

where X is distributed according to Fx(u) — y/u> with u € = [0, 1]. 

Now we discuss the L.L.N, for the current of first class particles that cross over the second 
class particle. Let £ be a configuration distributed according to ^00,0 and with a single second 
class particle at the origin - £ is represented on the left hand side of figure 2 below. Recall 
that the position of the second class particle at time t is denoted by X 2 {t). 

Denote by Jf r (i) the number of first class particles that cross over the second class particle, 
i.e. the number of first class particles at the right of Xa(t) at time t: 

Ji r (t)= £ (2.4) 

x>X 2 (t) 

As a consequence of Theorem l2.2[ we conclude that 

Corollary 2.3. Consider the tazrp starting from ^00,0 an d at the initial time put a second 
class particle at the origin and remove the first class particles there. Then 



J zr (t) 

lim — = (1 — Vx) 2 , in distributi 



and 



lim E u 

t-»+oo ^° 



where X is given in Theorem 



Now we consider a slight perturbation of the invariant state and we derive a L.L.N, for 
a single second class particle. 
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Theorem 2.4. Consider the tazrp starting from \i = u x + <5(_ co ]^ p -a with < A < p < oo 
and p — A < 1. At the initial time put a second class particle at the origin, let X 2 (t) denote 
its position at time t and let J 2 r {t) be the number of first class particles that cross over X 2 (t). 
Then 

lim — ^— = Z, in distribution under p*St 1 

t— >+oo t 

J zr (t) ( /~\ 2 * 
lim = 11 — vZ] , in distribution under u* St 

*->+oo t \ J 



and 



lim E^St 



Ji r (t) 



t 



= i + r log 



1 + A\ l 



A °\1 + P J (l + A)(l + p) 



where /!*(•) = p(-\X 2 (0) = 0), Z has distribution function given by Fz(u) = P ~"^ 1 A "^ , with 
m€3 = W{p), 0'(A)] and p{t,u) is the unique entropy solution of (2.1) with initial condition 
Pq X as in S2.3i . 

3. proof of Theorem |2.1| 

Proof. Fix a configuration ^ e N z and denote by J"(£) the current of particles that cross the 
time dependent bond ut during the time interval [0, t}. Then </"(£) is the number of particles 
of £ that are at left of the origin (including it) at time and are at the right of ut at time t, 
minus the number of particles of £ that are strictly at the right of the origin at time and 
are at left of ut (including it) at time t: 

t{x) 

Jt(0 = l {X^'((,)>ut} - HZ! («)<«*}' (3 - 1) 

i<0z=l x>Q z=l 

where for a site x, we label the t;(x) particles at this site (the first one being the one at the 
bottom and the t;(x) begin the one at the top) and we denote by the position at time t 

of a tagged particle initially at site x at position z. 

For x e Z denote by t x the space translation by x, so that for a local function / it holds 
that r x f(£) — /(r x £) and for a probability measure u e N z we have that J f{C)T x u(dti) = 
J f{ T x£,) n (dt;)- Now, we compute in two different ways: 



E[Jm}» P A<%) - J £[J t u (0]r-i/w(dO- 



Let £° and ^ be two copies of the tazrp starting from £° and £ x , respectively, such that ^° (5 1 ) 
is distributed according to p Pt x (t-iMp,a)- For any coupling p of p Pt \ and t_i^ p a and for any 
coupling P of $ and ^, last expression equals to 

dp(e,e)E[j?(e)-j?(t;% (3.2) 

where E is the expectation with respect to P. Obviously that in order to have a difference 
between the two fluxes, the configurations £° and must have at least a discrepancy at 
the origin at time 0, otherwise the difference is zero. Splitting this event by the number 
of discrepancies between the configurations and conditioning on D k which is the event that 
there are exactly k discrepancies at the origin between £° and £ x , last expression can be 
written as 



oo 

E 

fe=l 



dfi(z o ,e)E[j?(z°) - jne)\D k ]p(D k ). 



Now, we recall the dynamics of second class particles sharing the same site. Starting from 
D k , there are k discrepancies at the origin which we label from the bottom to the top as 
X 2 (0), ..,X k+ i(0). If the clock rings, the first particle to jump is X 2 (0) and then jumps X 3 (0) 
and so on. For j = 2, .., oo let Xj(t) denote the position at time t of the Xj(0) second class 
particle. Since there are k discrepancies, the difference between the currents can vary from 
1 to k, depending on the relative positions of these second class particles with respect to the 
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moving bond ut. If j is the difference between the fluxes, then we must have Xj +1 (t) > ut but 
Xj+ 2 {t) < ut. Notice that having Xj + i(t) > ut, it implies that X 2 (t) > ut, ..,Xj(t) > ut, since 
the Xj + i(t) -second class particle only jumps to the right if the particles with less degree of 
class are at its right. This is the key point in the proof where we use the total asymmetry of 
jumps. In the presence of partial asymmetry the second class particles do not preserve their 
order. This fact is crucial for our conclusions and allow us to write last expression as: 



oo k— 1 
fe=l j=l 



K,— ± 

Y,3P^{X j+1 {t) > ut,X 3+2 (t) < ut) + fcP fe (X fe+1 (t) > ut)) (ji(D+) - HD-)), (3.3) 



where fi k is the coupling measure conditioned on D k and D% = £ x ) : £°(0) = ^(0) + 
(— )k}. In D^(D^) there is a contribution for the current with a positive (negative) sign, since 
the current can vary from 1 (—1) to k) as long as the discrepancies at time t are at the 
right of ut Since we are dealing with Geometric product measures, it is not hard to show that 



m k 



P(D+) - M(D fe ") = (: 



1 + P 



A 



1 + A 



(3.4) 



.l + p + A> 

and in sake of completeness we prove in Remark l3.ll In order to keep notation simple, write 
down the probability {X. J+ i{t) > ut, X J+2 (t) < ut) as pk(j, t). Then H3.3D can be written as 



fc=i 



oo fc— 1 

J2 CE,jPkU,t) + kP Pk {X k+1 {t) > ut))m fe . 

3=1 



(3.5) 



Now, notice that 



fe-i 

E 

3 = 1 



fe-1 

jp k (j,t) = (P^(X J+ i(t) > ut) - (fc - \)P fik (X k+1 (t) > ut)), 



3 = 1 



since by a simple computation p k (j, t) can be written as P, Ik (Xj + i(t) > ut) — P IIk (Xj +2 (t) > ut). 

Last equality follows from the fact that if {Xj +2 (t) > ut} then we must have for sure 
that {X 3+ i(t) > ut}. Again, total asymmetry of jumps in invoked. If the jumps are partially 
asymmetric and the difference of currents is j, then in D k j of the k discrepancies are at 
the right of ut. This can happen a number of Cj 5 events but we do not have control over the 
probabilities of each one of them, and for that we cannot get to a similar result. Collecting 
these facts together, ( 13.51 ) can be written as 



Y / m k ( K Y / P^(X J+ i(t) > ut) 



k=l 



3 = 1 



Notice that Pfi k (Xj + i(t) > ut) does not depend on k, since this is the probability that the (j + 
l)-th second class particle at the origin (with other k second class particles at the same site), 
having crossed the line ut at time t and that does only depend on j and t. This observation 
allows us to apply Fubini to change the order of summation in last expression to write it as: 



By (US it holds that 



So far we have that: 



J2( p fi( X 3+i(t)>ut)J2mk 

3 = 1 k=j 



J2 mk = T 



k=3 



+ P + a 1.(1 + py- 1 (i + xy- 



dn(e,e)E[jne) - jne)} = £ (^pw*) > 



i 



3 = 1 



+ P+X 



(3.6) 



/V 



(1 + p)^- 1 {i + xy- 
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Now, compute A3. 2D by coupling p pX and T_i/i Pt x in such a way that £ x = r_!^°. By the 
definition of the current (see d3.ll >) we write J"(^°) — -/"(C 1 ) as 

C(x) t°(x) e(x) 

X X 1 {^ x ' i (e o )>«t} ~zJ X 1 w^«°)<«*} "EE 1 {*r j, (€ i )>«t} + zZ X ^-'(ei)^}' 

K<0 2 = 1 £<0 2 = 1 X<0 2 = 1 X<0 Z—l 

where and Xf' z (^ 1 ) denotes the position of a tagged particle initially at site x at 

position z for the process starting from £° and £ x , respectively Since i; 1 = t_i£° and writing 
the expression above in terms of the process £° we get to 

X X ~ 

xez 2=1 

Recall that £° is initially distributed according to ^ p ,a- Applying expectation with respect to 
p P: x to last expression, we obtain that ( 13.21 1 equals to 

E, p J$(ut + l)}-E, p je(l)]. 

Since £°(1) is distributed according to p\ it follows that E tlp x [£°(1)] = A. On the other hand, 
by the convergence to local equilibrium (see 13]) it follows that: 

lim E^ x [$(ut+l)] = p(l,u), 

where p(t, u) is the entropy solution of d2. ID with initial condition pft as in d2.3D . Putting 
together the previous computations the proof ends. □ 

Remark 3.1. Here we compute p(D^), i.e. the p probability of having k discrepancies between 
£° and such that £°(0) = ^(0) + k. Recall that £° and S, 1 are distributed according to p Pt \ 
and T-ipp^\, respectively. Notice that 



P(Dt) = fi( (J {e°(0) = n + k,e(0) = n}) 

CO 

E/ p y i + k 1 
VTT 



n=0 



1 + /A1 + AV 1 + A 
1 ( P \ k 



(1 + p + A) \l + pj 

In the second equality above, we use the fact that p is the coupling measure and it is a product 
measure with marginals p Pi \ and T-ip Pt \. A similar computation shows that 

1 /A 



(1+p + X) VI + A 

Remark 3.2. Last result also holds for totally asymmetric jumps given by <?(•) as in the 
beginning of section \2.1\ This is a consequence of the attractiveness property and the local 
equilibrium convergence for the process starting from p p ,\, see (3J. 



4. TAZRP STARTING FROM /i^o 

Consider the tazrp starting from poofi, put a second class particle at the origin and remove 
the first class particles there, i.e. the process starts from £ introduced before the statement 
of Corollary 12.31 We prove a L.L.N, for the position of the second class particle and for the 
current of first class particles that cross over the second class particle. We start by the former. 
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4.1. Proof of Theorem \2M 

Proof. Let | be the configuration with infinite particles at negative sites, the origin and all 
positive sites empty: 

oo, if z < —1 
0, ifz>0 



Recall the definition of J"(£) from the proof of Theorem |2. li For £, it follows that 

X>1 

i.e. J"(£) is the number of first class particles at the right of ut at time t. By the Kolmogorov's 
backwards equation we have that 



(4.1) 



Notice that since at site —1 there are infinite particles, there is only one discrepancy at 
the origin between and £. In the figure below the configuration represented on the top 
is £ and the other one is £ -10 . 



ooooooo 



FIGURE 1. One single discrepancy at the origin. 



Now we consider to copies £ t ° and £ t x of the tazrp starting from £ -10 and £, respectively, 
and we use the basic coupling on them. The basic coupling was introduced in Ifl2ll and it 
means that, by attaching a Poisson clock of parameter one to each site of Z, when the clock 
rings the particle jumps to the right neighboring site. With this coupling both processes 
use the same realizations of the clocks. By the attractiveness property together with the 
conservation of the number of particles, there is still a unique discrepancy between £ 4 ° and £ t x 
at any time t > 0. Let X 2 (t) denote the position at time t of this discrepancy and notice that 
by construction X 2 (0) = 0. Now, since J"(£~ 1,0 ) — = l{x 2 (t)>«t}> then d4. It is equal to 

P(X 2 (t) > ut). 

On the other hand, by writing J t n (£) as a martingale plus a compensator: 



J t "(0 = M t (£) 



(«*)>!} dS ' 



we get that 



At^oo.o [Jt(Q] - En a 



1, 



«t(«t)>i}J' 

since the mean of martingales is constant and the martingale M t vanishes at time 0. Now, 
by the convergence to local equilibrium it follows that 



lim En 

t->+oo 



1 {6(«*)>i} 



= </>(p(l,u)) = 1- 



Here p(t, u) is the unique entropy solution of d2.ll * starting from pfr x as in 02.3D with p = oo 
and A = 0. Collecting these facts together the proof ends. □ 
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Now, we prove the L.L.N, for the current of first class particles that cross over the second 
class particle. 

4.2. Proof of Corollary |2~3l 

Proof. Recall from Theorem 12.21 that X2{t)/t converges in distribution as t — > +oo to X dis- 
tributed according to Fx{u) = ^/u, for u e 3x- By the definition of Jf r (£), see ( 12.41 ); from the 
hydrodynamic limit (see section [272] > and Theorem l2.2[ it follows that 

lim 2 y ' = / p(l,u)du= / (-=-l)du= (1- Vx) , 

t^oo t J x J x \y/U ) \ J 

where p{t, u) is the solution at time t of ( 12.11 1 with initial condition pg' A as in 02.31 ) with p = oo 
and A = 0. This finishes the proof of the first claim. For the second, since X has distribution 
given by F x (u) = yfu for u e [0, 1], then 

lim E,, 



Mo 

t-V + OO ^ 



\Ji r (ty 


-I 


i t - 





» 1 -.du= 1 - 



3' 

□ 



5. Law of Large Numbers more general p 
5.1. Proof of Theorem i3 

Proof. The proof follows the same arguments as in the proof of Theorem 1 of (H. The main 
features of the process that we need in order to achieve the result are: the attractiveness 
property, the existence of at most one discrepancy when applying any coupling to p and r_i/i 
and the convergence to local equilibrium for the process starting from p. 

As in the proof of Theorem 12.11 fix £ e N z and denote by </"(£) the current of first class 
particles that cross the moving bond ut. We compute in two different ways: 

J E[J?(t)Md£) - J E{J?(0}r-ip(dO- 

Let £° and £ t x be two copies of the tazrp starting from £° and £ x , respectively, such that £° (£ 1 ) 
is distributed according to p (r-xp). For any coupling p of p and T-\p and any coupling P of 
£ f ° and last expression equals to: 

where £ is the expectation with respect to P. The measure p induces in p the property 
that £°(a:) = fo r i^O and at the origin £° has an extra particle relatively to f 1 , with 
probability p — A. Let X 2 (t) denote the position at time time of this discrepancy. Since there 
is at most a discrepancy between £° and the difference between the currents is at most 
one, and that happens as long as X 2 (t) is at the right of ut at time t. With this observation 
d5.ll ) can be written as 

(p-X)P^s t (X 2 (t) >ut), 

where p* S t is the distribution of the process at time t starting from p* . 

On the other hand, we couple p and T-ip in such a way that = t_i£°, which allows to 
rewrite ( 15.1D as 

E^(ut+l)]-EX0)]- 
Since under p, is distributed according to p\, then £^[£°(1)] = A. Now, by the conver- 
gence to local equilibrium (see remark [§31 ), it follows that: 

lim E^iut+l) = l] = p(l,«), 

where p(l, u) is the unique entropy solution of ( 12.1D with initial condition Pq' a . This is enough 
to conclude the first claim of the Theorem. 
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Now we derive the L.L.N, for the current of first class particles that jump over the second 
class particle. Denote by Xi (t) the position at time t of the first class particle initially at site 
— 1. Since first class particles do preserve their order, the current of first class particles that 
cross over the second class particle from time to time t, is equal to the number of first class 
particles at time t between the sites Xi(t) and X 2 (t). Since X%(t) is the first to jump over the 
second class particle, then 

x 1 (t) 

Ji r (t)= E 

x>X 2 (t) 

Notice that at positive sites, is distributed according to So that a L.L.N, for X\{t) 
holds (see IflSH): 

lim « = 1 A = 1 
t^+oo t 1 + A 1 + A 

On the other hand, from the hydrodynamic limit and the previous result we get that: 

lim 2 - - = / p(l,u)du = / ^—du+ f Xdu. 

t Jz Jz J^-z 

Now, a simple computation ends the proof. For the second claim, it is enough to notice that 

r J zr (i\ i r<t>'( x ) 1 

lim E,, St \ J -lM\ = (1 - -r^du. 

□ 

Remark 5.1. We recall that the hydrodynamic limit for the tazrp and for more general 
jumps rates given by g(-) (as defined in section \2J) starting from a, can be obtained by the 
Relative Entropy method following the same arguments as done for the tasep in ©. As a 
straightforward consequence we obtain the local equilibrium convergence. For details see 
chapter 6 of [10]. In possession of last result, we can obtain for the zero-range with jumps 
rates given by <?(•) the same result as in Theorem \2.4\ but with Z is distributed according to 

Fz(u) = p -%™ r \ for u e (g'(A))- 1 ]. 

Remark 5.2. We remark here that last result also holds for partially asymmetric zero-range 
processes with jump rate p(-) with positive and finite mean: < 7 := J2 x ei,P( x ) < 00 ■ The 
proof follows the same lines as in the proof of Theorem \2.4\ The main difference is that p(t, u) 
is the solution of the corresponding hydrodynamic equation. 

6. Coupling with tasep 

6.1. Law of Large Numbers. In this section we reprove Theorem l2.2l and Corollary I2.3l bv 
coupling the tazrp with the tasep. Recall that we consider the tazrp starting from /ioo,o, with 
a second class particle at the origin and no first class particles there. We start by showing 
a L.L.N, for the position of the second class particle initially at the origin. First we recall a 
result that we will use in the sequel: 

Theorem 6.1. (M) 

Consider the tasep starting from v Pi \ with < A < p < 1. At time zero put a second class 
particle at the origin regardless the value of the configuration at this point and let Y 2 (t) denote 
its position at time t. Then 

lim 2 ^ - = U, in distribution 

t— >+oo t 

where U is uniformly distributed on [(1 — 2p), (1 — 2A)]. 

The proof of last result was given in (5]] and extended for more general asymmetric exclu- 
sion processes in |4|. By coupling the tazrp with the tasep, we can show a L.L.N, for X 2 (t) 
and show that: 
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Proposition 6.2. Consider the tazrp starting from and at time zero put a second class 
particle at the origin and remove the first class particles there. Let X 2 (t) denote the position 
at time t of this second class particle. Then 

Urn — ^— = f — — \ , in distribution (6.1) 



t— >+oo t 

where U is uniformly distributed on [— 1, 1]. 

Proof. The idea of the proof is to couple the tazrp £ t with the tasep n t and to establish a 
relation between the position of the second class particle on the tazrp with some microscopic 
function on the tasep. Now we explain the coupling we use between the two processes. 

Suppose to start the tazrp from £ - the configuration distributed according to and 
with a single second class particle at the origin, see the left hand side of figure 2. On the 
other hand, start the tasep from uifi, remove the first class particle at the origin if necessary 
and add there a second class particle. On the tasep, initially label the first class particles by 
denoting the position of the first class particle at site — i at time by Xj(0). Recall that the 
occupation variables in the tazrp are represented by (f(a:))xez- We relate both processes in 
such a way that in the tasep the distance between two consecutive first class particles minus 
one, becomes the number of particles at a site in the tazrp representation. With this in mind, 
we define for i > 1: £ (i) = x i+ i(0) — 2^(0) — 1, at the origin both have a second class particle 
and for i < — 1, £ (i) = oo since on the tasep there are no other first class particles at the 
right of cci(0). On the left hand side of the figure below, we represent the tazrp, the second 
class particle is represented by ®, while first class particles are represented by • and holes 
by o. On the right hand side, the tasep is represented, first class particles are represented by 
• and holes by o. 



ISOOOOOOO •••••••^f^OOOGQQO 



FIGURE 2. Coupling between tazrp and tasep 

At first notice that in the tazrp, a first class particle can jump to the site where the second 
class particle stands and the latter does not move. This corresponds to a first class particle 
jumping over the pair o« on the tasep representation. On the other hand, in the tazrp when 
the second class particle jumps to the right, this corresponds to a jump of the pair o« over 
a hole on the tasep. With this relation, the position of second class particle at time t in the 
tazrp, corresponds to the number of holes to the left of the pair o« at time t in the tasep. 
Now, we also notice that for the tasep, the dynamics of the pair o» in the configuration above 
is the same as the dynamics of a second class particle at the origin in the configuration 
with all negative sites occupied and all positive sites empty, see |6] for details on this last 
correspondence. 

In order to summarize the established relations we introduce some notation. Recall that 
X 2 (t) denotes the position of the second class particle at time t for the tazrp starting from £. 
Let Jf e (£) (Jf(i)) denote the number of first class particles at the right of the second class 
particle at time t in the tasep (tazrp) and let fff e W denote the number of holes to the left of 
the second class particle in the tasep at time t. For this particular initial conditions, by the 
relations mentioned above we have that 

Ji e (t) = Jf(t) and H^{t)=X 2 {t). (6.2) 

So, the position of the second class particle in the tazrp starting from £, corresponds to the 
number of holes at the left of the second class particle for the tasep starting from v lfi with a 
second class particle at the origin. Recall the L.L.N, for H^it) shown in |7], namely that 

lim — 2 - - = ( — — — 1 , in distribution (6.3) 
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where XL is a Uniform random variable on [—1,1]. Then, the equality on the right hand side 
of ( 16. 2D together with last result finishes the proof. So, we identify X in the statement of 

Theorem EJ as {^r) ■ □ 

Remark 6.1. We remark that the convergence in H6.3i in fact takes place almost surely. As a 
consequence, the convergence of the second class particle in S6.1i also holds almost surely. 

Now we reprove Corollary |2.3l bv identifying the limit random variable for which J£ e (t)/t 
converges. 

Corollary 6.3. Consider the tazrp starting from and add at time a second class parti- 
cle at the origin and remove the first class particles there. Then 

J| r (i) /l - U\2 
lim = — - — , almost surely 

4-S-+0O t \ 2 / 

where XL is uniformly distributed on [—1, 1]. 

Proof. By d6.2| > and since it was shown in (7) that J| e (i)/i converges almost surely to (ht^) 
where XL is uniformly distributed on [—1,1], the result follows. □ 

6.2. Crossing probabilities for a second and a third class particle on the tazrp. 

Corollary 6.4. Consider the tazrp starting from /ioo,o ond at time zero add a second class 
particle at site and a third class particle at site 1 and remove the first class particles there. 
Denote by X 2 (t) (X 3 (t)) the position at time t of the second (third) class particle. Then 

lim P(X 2 {t)>Xz(t)) = \. 

Proof. The proof of this result is a consequence of the relation between the tazrp and the 
tasep that we introduced above and of Theorem 2.3 of (4). The probability that we want 
to compute is the same as starting the tazrp from a configuration with all negative sites 
with infinite particles, a single second class particle at the origin and another single one at 
site 1 and all positive sites empty (see the left hand side of the figure below), and ask for 
the probability of the second class particle at the origin to jump over the other second class 
particle. In the picture bellow, on the left hand side the tazrp is represented and first class 
particles are denoted by •, the second class particles by ® and holes by o. On the right hand 
side the tasep is represented and first class particles are represented by • and holes by o. 



Figure 3. Relation between tazrp and tasep with a second and a third class particle 



As explained in the proof of Proposition 16.21 the dynamics of a second class particle in the 
tazrp corresponds to the dynamics of the pair o» in the tasep. In this case, the dynamics of 
the two second class particles in the tazrp representation corresponds to the dynamics of the 
two pairs o»o» in the tasep. But, the dynamics of the two pairs in the tasep, is the same as 
the dynamics of two second class particles for the tasep starting from a configuration with 
all negative sites occupied, two second class particles (one at the origin and other at site 1) 
and all positive sites empty, see [|6]| and Hi for details. Also, notice that this last dynamics 
corresponds to the dynamics of second and third class particles for the tasep starting from 
the configuration with all negative sites occupied, all positive sites empty, a second class 
particle at the origin an a third class particle at site 1, see (4). As a consequence, the claim 
above follows straightforwardly from the result in Theorem 2.3 of (4). O 
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